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Spatially inhomogeneous shear flow occurs in entangled polymer solutions, both as steady state shear banding 
and transiently after a large step strain or during start up to a steady uniform shear rate. Theoretically, steady 
state shear banding is a hallmark of models with a non-monotonic constitutive relation between total shear 
stress and applied shear rate, but transient banding is sometimes seen in fluids that do not shear band at 
steady state. We model this behavior using the diffusive Rolie-Poly model in a Newtonian solvent, whose 
steady state constitutive behavior can be monotonic or non-monotonic depending on the degree of convective 
constraint release (CCR). We study monotonic steady state constitutive behaviour. Linear stability analysis 
of start up to a sufficiently high shear rate shows that spatial fluctuations arc unstable at early times. There is 
a strong correlation between this instability and the negative slope of the (time dependent) constitutive curve. 
If the time integral of the most unstable eigenvalue is sufficiently large then the system exhibits transient shear 
bands that later vanish in steady state. We show how perturbations, due to fluctuations or the inhomogeneous 
stresses, can trigger this instability. This transient behavior is similar to recent observations in entangled 
polymer solutions. 



I. INTRODUCTION 



ticularly in the response of highly entangle d poly meric 



fluids to a step strain . Osaki and Kuratal ( 198(1 ) and 



The original theory of iDoi and" Edwardsl (|l989t) pre- 
dicts shear banding in entangled polymers, when the 
shear rate exceeds the reciprocal of the timescale Td 
for one-dimensional diffusion (reptation) within an effec- 
tive mean field tube of constraints from the surrounding 
molecules. For shear rates 7 > 1/Vd, DE theory pre- 
dicts a decreasing shear stress a xy (j) due to enhanced 
tube alignment, and thus leads to a non-monotonic con- 
stitutive curve <Jxy(j) an d an instability to an inho- 
mogeneous state in which two states, or 'shear bands', 
of different viscosi t y and molecular alignment coexist 
( Mc Leish and Ball . Il986|) . However, banding was not 
inferred at that time in polyme rs; for example, the ap- 
paren t viscosity measurements of lMenezes and Graesslevl 
( 19821 ) on entangled polymers are consistent with a 
weakly increasing shear stress for jTd > 1, rather 
than the stress plateau that is expected from band- 
ing. By contrast, banding is preval en t in wormlike 
micelles dRehage and Hoffmann!. 1 199 it iBerretl. 120051: 
ICates and Fielding). 120061: lOlmstedl . 120081) . This was ex- 
plained by ICated (|199oT k who combined reptation with 
micellar breaking to predict shear banding at stresses 
and shear rates in good qualit ative agreement with ex- 
periments ( Spenlev et oUll993) . It is now realized that a 
strictly constant stress plateau is not obtained for band- 
ing in a curved geometry, such as cone and p late or 
cylindrica l Couette fixtures (lOlmsted et aZ.l . r2000h . so the 
results of enezes and Graesslevl (|l982r )mav in fact be 
consistent with shear banding. 

Discrepancies with DE th eory have exis t ed vir tually 
since the theory's inception ( Fukuda et all , Il975f ). par- 



IVrentas and Graesslevl ( 1982f ) observed stress relaxation 



that was much faster than predicted by DE theory. This 
is particularly prevalent for entanglement numbers Z > 
60 for relaxation after large step strain ()Osaki Il993[) . 
and ha s been t ermed anomalo us or "type C" relaxation 
" Osaki Il993fc IVenerud . I2005D . iMarrucci and Grizzutil 
19831 ) (MG) pointed out that the DE theory contains an 



elastic instability; they showed that the free energy can 
be a convex function of strain for large enough step strain 
(for jTd 3> 1). This would then lead to an elastic insta- 
bility and thus inhomogeneous (shear banding) behavior. 
They showed how this in stability can persis t even when 
one includes relaxation. iMcLeish and Ba ll ( 1986h later 
suggested that the DE theo ry could explain the measure- 
ments of the spurt effect bv lVinograd"ovl(ll973l). in terms 
of the inherent instability to shear banding. Kolkka et al\ 
( 19881 ) used the Johnson-Segalman model to study the 
nature of mechanical instabilities inherent in fluid with a 
non-monotonic constitutive relation. It was later demon- 
strat e d that wall sl ip plays a major role in spurt ( Wand . 
119991 ; iDennl . l200l[ ). iMorrison and Larson! (| 19921) tested 
Marrucci and Grizzuti's idea by using entangled solu- 
tions at different concentrations to mimic the effect of 
relaxation. They concluded that the Marrucci- Greco 
mechanism gives a qualitatively correct description of the 
anomalous relaxation, which is consistent with incipient 
inhomogeneities. More r ecently, Venerus ( Venerud . 120051 : 
IVenerus and Nairl . [2006^) studied step strain experiments, 
and concluded that many of the anomalous "type C" re- 
sponses could be explained by experimental conditions 
such as wall slip or rheometer compliance; however, a 
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number of experiments could not be accounted for in this 
way. 

In the last decade or so the advent of high resolu- 
tion velocimetry has led to a series of exciting experi- 
ments, principally by Wang and co-workers, that demon- 
strate evidence of shear banding i n entangled polymers 



dTapadia and Wand. 120031 12004: IBoukanv and Wan" 



l2007t iHu et all l2008at IRavindranath et all |2~ 



Very well-entangled solutions, with entanglement num- 
ber Z > 40 — 50, displa y well- defi ned stead y state 



shear bands (IWang et cd 



2008a 



IRavindranath et al 



Hu et all 12007 



IBoukanv and Wang 



2009bJ). This occurs in both synthetic poly mers, 



such as polybutadiene (iRayindranath et~all I2008D. and 
DNA (IBoukanv et allhooA IBoukanv and Wand . l2009at 



IHu et aUl2008bl ; IBoukanv and Wand . l2009bl) . Polvdis- 
perse systems e xhibit smoother, less well-d efined (or no) 
shear banding ( Boukanv and Wanel . 120071 ). possibly be- 
cause the ma ny timescales smear out the instability (as 
suggested by iDoi and Edwards! (|l979h in their original 
paper). 

There have been numerous improvements upon the 
original DE theory, nota bly to incorpora t e convective 
const r aint release (CCR) ( Marruccl Il996t iMead et~al\ , 
119981 ; iMilner et ai] . l200ll ). CCR can reduce the sever- 
ity of t he DE instabili t y and potentially render the fluid 
stable (jGraham et all \200S\ ), and DE-CCR models can 
capture, at least qualitatively, many of the shear banding 
sig natures of these exper i ments. This was recently shown 
by lAdams and Olmsted! ( 2009bllah in calculations base d 
on the Rolie-Poly model ( Likhtman and Grahaml I2003T) , 
a simple one-mode differential version of DE theory that 
incorporates CCR and chain stretch. However, there is 
no consensus yet as to the correct level of CCR required 
to describe existing experiments. 

In related work, IZhou et al. (2008) showed that a non- 
monotonic constitutive relation base d on the "par tially 
extending convective strain" model of lLarso n1 (|l984l) also 
reproduces many of the features of the recent experi- 
ments, including recoil during startup onto the banding 
plateau, and strain localization immediately after a step 
strain effected by shearing at a very fast finite shear rate. 

Numerous experiments also show inhomogeneous be- 
havior in fluids that do not display steady state shear 
banding (smaller entanglement number Z). Examples 
include: (1) transient band formation and recoil dur- 
ing startup, which eventually gives wa y to homog e neous 
shear flow dTapadia and Wand. 120061: IHu et all l2007t 



IRavindranath et all 120081 : IBoukanv and Wand . l2009bl) : 
(2) extremely sharp shear banding during portions of 
the cycl e during Large Amplitude Oscillatory Shear 
(LAO S ) dTapadia et atl.l2006l : IRavindranath and Wanel . 
120081 : IZhou et all l2010f) : (3) inhomo geneous response 
(including negative velocity recoil ) after step strain 

performed at high she ar rates (|Wang et all 120061: 

IRavindranath and Wanel 120071) . 

At a constitutive level, one cause of steady state shear 
banding is a non-monotonic steady state constitutive 



relation (such as the DE theory with sufficiently weak 
CCR), which cannot support stable homogeneous steady 
states for a range of shear rates. However, the same cri- 
terion clearly cannot determine the apparent instability 
to transient banding in fluids that don't shear band in 
steady state (i.e. that have monotonic flow curves). Evi- 
dently these fluids are dynamically unstable, which need 
not be related to a steady state non-monotonic constitu- 
tive relation. 

Indeed, pronounced t ransient shear band ing has re- 
cently been reported by iDivoux et al\ (|2010l ) in a yield 
stress fluid with a monotonic steady state constitutive 
curve, and captured theoretically in shear transfo rmation 
zone (STZ) theories (|Manning et all 120071120091^ . a mod 



ified soft glassy rheology (SGR) model (Moorcrof t et al 
2011 ). a simplified fluidity model dMoorcroft et al 



2011 ). and a mesoscopic model of plasticity ( Jaglal 20101) 



Several, possibly related, explanations for the transient 
banding can be envisioned: 

1. The dynamical equations of motion passes through 
a regime of parameter space, as a function of time, 
in which the homogeneous instantaneous state is 
dynamically unstable to small spatial perturba- 
tions. Such perturbations will initially grow in 
time before eventually decaying to a homogeneous 
steady state. However, if the instability is strong 
enough the perturbations could grow into a macro- 
scopically observable transient band before decay- 
ing. We will perform this linear stability calculation 
here, and show that this gives an understanding 
of the transient homogeneities. A simil a r calc ula- 
tion was done bv lFielding and Olmsted! ( 20031 ) for 
wormlikc micelles undergoing steady state banding. 

2. The fluid possesses, at any observation time t m , an 
instantaneous constitutive curve a xy (7, t m ). Such 
a constitutive curve can be constructed by perform- 
ing a number of startup evolutions (calculations or 
experiments) at different shear rates, which then 
define the loci of stress as a function of shear rate 
for given observation times t m after flow inception. 
Since most of these points are not in steady state, 
the resulting curve need not be monotonic even if 
the shear rate remains homogeneous. Then, one 
could associate a negative slope daxy^^t^/dj < 
with transient inhomogeneities. However, this 
only addresses a subspacc of the full parameter 
space considered by the first scenario above and 
could thus only serve as a rough guideline. One 
must also carefully distinguish between theoretical 
calculations in which a homogeneous shear rate can 
be specified, and the experimental protocol above, 
in which the fluid will be expected to become in- 
homogeneous after the onset of an inst ability. This 
proced ure was recently suggested by lHaves et al\ 
( 20101) in experiments using parallel plate rheome- 
try of polymer solutions. 
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For very strong shear rates these viscoclastic flu- 
ids respond elastically, like a non-linear solid, and 
can display a stress overshoot. In t his ca se the 
mechanism of iMarrucci and Grizzutil (Il983h might 
be expected to apply, as follows. For a given im- 
posed shear rate 7 the stress a xy is a function of 
strain 7 = jt. A solid with a negative gradient 
da xy /dj < has an effective negative differen- 
tial shear modulus, which leads to clastic instabil- 
ity. This would imply that, for shear rates large 
enough to remain in the elastic regime, a stress 
overshoot upon startup could signify an instabil- 
ity to inhomogeneous flows. [Note that the MG 
argument is strictly for an ideal step strain, rather 
than a strain incurred during a fast but finite shear 
rate] . ISui and McKeniial ( 20071 ) pointed this out in 
their recent vi sualization study of the sa me materi- 
als studied bv lTapadia and Wand (120061) . This sce- 
nario is suggested by data showing correlations be- 
tween transient bandi ng and stress overshoot s (e.g. 
the Figures 3 and 4 of iBoukanv et all l|2008j) ). and 
by calculati ons that show similar correla tions (e.g. 
Figure 3 of I Adams and Olmsted! (j2009bl) . for both 
monotonic and non-monotonic constitutive curves). 



In this paper we explore these different scenarios, and 
thus the detailed conditions necessary for transient inho- 
mogcncities, using th e diffusive Rolie-Poly (DRP, or RP) 
model as an example ( Likhtman and Grahaml . l2003|) . We 
calculate its instantaneous linear stability for parameters 
ranging from monotonic to non-monotonic constitutive 
behavior. We show that even fluids with monotonic con- 
stitutive curves, which have stable homogeneous steady 
states, can have periods of instability during startup flow. 
For some parameters, strong linear instability is shown 
to lead to transient banding and negative velocity recoil 
in the full non-linear dynamics. 

Since instabilities require an initial perturbation to 
manifest in eventual non-linear growth, we will study 
two natural perturbations: (1) stress gradients, as found 
in typical circular rheometric devices (cone and plate or 
cylindrical Couette), and (2) non- uniform (noisy) initial 
conditions. We find a strong link between transient in- 
homogeneities, stress overshoots, and an instantaneous 
non-monotonic constitutive curve, but we leave further 
detailed analysis of this for future work. 

We compare our calculations with recent experiments 
on transient data after startup to steady state. For con- 
ciseness, we do not study the inhomogeneous response 
to a step strain, though our results give a quantitative 
method for understanding this behavior as well. 



II. THE DIFFUSIVE ROLIE-POLY MODEL 



A. Momentum Balance 



Newton's force balance for the fluid is given by 



Dt 



(1) 



where p is the density, a is the total stress, and D/Dt is 
the material derivative. We will use the full equation set 
(at finite but small Reynolds number) for the linear sta- 
bility analysis in section lTV A[ but take the zero Reynolds 
number (creeping flow) limit for calculating the full spa- 
tially resolved non-linear dynamics in section IIVBI In 
the latter case the equation of motion reduces to 



V • cr = 0. 



(2) 



For the experimental cases of interest we expect very 
small Reynolds number Re <C 10~ 3 . 



B. Stress 

The total stress cr in a polymer solution is assumed 
to comprise an elastic stress carried by the backbone of 
the polymers (due to their stretching and orientation), 
and viscous drag against solvent and other polymers. In 
polymer melts the viscous stress ca n usually be safely 
negle cted, particularly for weak flows (|Doi and Edwards! 
1989) . However, for strong flows, where DE theory pre- 
dicts a strongly decreasing polymeric stress due to tube 
alignment, it is necessary to incorporate additional vis- 
cous stresses from the faster degrees of freedom. This 
includes the Newtonian solvent viscosity of solutions, as 
well as fast Rouse modes and intcrpolymer viscous fric- 
tion neglected in the simplest tube models. In addition 
to this crucial argument, a second contribution to the to- 
tal stress is physically necessary to describe steady state 
shear banding in planar Couette flow, in which the shear 
rate is inhomogeneous but the total stress must be ho- 
mogeneous. 

Hence, we represent the total stress as two separate 
components: fast Newtonian (or solvent) degrees of free- 
dom, and the slow viscoelastic stress GW: 



cr = -pi + 2r?D + GW, 



(3) 



where I is the identity tensor, D = ^ [Vv + (Vv) T ] , v is 
the velocity field, p is the isotropic pressure determined 
by incomprcssibility (V-v = 0), 77 is the solvent viscosity, 
and G is the plateau modulus. In this representation 
the quantity W is the polymer, or viscoclastic, strain, 
whose stress GW is parametrized by the elastic modulus 
G. Together with a DE-like constitutive relation for W, 
Eq. can yield either a monotonic or a non-monotonic 
constitutive relation, with a Newtonian high shear rate 
branch. 
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FIG. 1. Constitutive curves for the Rolie-Poly model, show- 
ing monotonic and non-monotonic behavior as a function of 
entanglement number Z\ for j3 = 0.65 and e — 10 . 



C. Rolie-Poly Model 

There are many possible constitutive mod e ls for 
the polym e ric st rain W ( Islam and Archeil 1 2 ll : 
iMead et all Il998l) . but here we will use the diffu- 
sive R olie-Poly (DRP) model of iLikhtman and Grahaml 
(2003). This is a single mode a pprox i mation to the 
GLAMM model (iLikhtman et all boOOl : iMilner et all , 
l200ll ; iGraham et aZ.l . l2003 ). and includes convective con- 
straint release (CCR), reptation of the polymers within 
their tubes, and the stretching of the polymer chains. 
The constitutive equation for the deviatoric part of the 
viscoelastic strain W is 



(d t + v • V)W - (Vv) T W W (Vv) + — W = 

Td 

2D - — (1 - A) (I + W + (3A- 25 W) + £>V 2 W, (4) 



where 



A(W) = (1 + trW/3) 



-1/2 



(5) 



Td is the disengagement time, tr is the Rouse time, and /3 
parametrizes the efficiency and rate of CCR. The ratio of 
disengag ement to Rouse tim es defines the entanglement 



number ([Larson et ~al\. I2003D 



3 t r 



(0) 



We have added 'diffusion' (the term with coefficient V) 
to the original Rolie-Poly model, to be able to resolve 
spatial structure dur ing shear banding ( Lu et all l2000t 
lOlmsted et al l l200fl The width of the interface between 
shear bands is proportional to y/V. 

The RP model does not reproduce all the same consti- 
tutive behavior features of the full GLAMM constitutive 
model (notably, the high shear rate behaviors of the nor- 
mal and shear stresses differ), primarily because only a 



single mode of the GLAMM model has been kept in ob- 
taining the RP model. However, it is simple enough to 
perform spatially resolved simulations in the presence of 
the added diffusive term. The additional solvent stress, 
277D in Eq. ([3]), was not included in the original formula- 
tion, which was focused on melts and not on shear band- 
ing behavior. 



D. Choice of Parameters 



Likht man and Grahaml (|2003h chose /3 = 0.5 and 5 = 
— 1/2, which optimized the comparison with transient 
and steady state measurements within the assumption 
of homogeneous flows. The negative value for 5 ensures 
that CCR decreases for large stretch, due to the increased 
number of effective entanglements. The high shear rate 
scaling of the total shear stress a xy in the RP model is 



1 



G 



' 3 

s 6Z/32 A 



(6<0) 



+ ejT d (<S > 0), 



where 



_n_ 

Gr d 



(7) 



(8) 



is the ratio between the solvent viscosity and that of the 
quiescent entangled solution. Typical experimental val- 
ues f or the parameters are e ~ 10~ 3 — 10~ 5 and Z ~ 15 — 
200 (lAdams and Olmstedl . l2009bl iTapadia and Wand . 
120061) . 

As can be seen above, the choice S = —1/2 leads to a 
Newtonian high shear rate branch, with slope dominated 
by the entanglement contribution for typical values of 
6 and Z. This is incons istent with the experiments of 
Tapad ia and Wand ()2003D . which suggest a scaling a xy ~ 
7 1 / 2 at high shear rate. This scaling can be produced by 
choosing (5 = 1/2 (Eq. [7]), which leads to 



u xy 

~G 



3^ 



ITd 



6Z/32 



1/2 



£7T d . 



(9) 



This has an intermediate scaling a xy ~ 7 ' , before giv- 
ing way at the highest shear rates to Newtonian behavior. 
However, this choice does a poorer job at th e valid ation 
tests carried out by ILikhtman and Gr ahaml (l2003h and 
has the incorrect physical interpretation of the effect of 
stretch on CCR (as noted above), and for these reasons 
we use 5 = — 1/2 here. 

The Newtonian power law in the high shear rate 
branch, which applies for our choice S = —1/2, produces 
a relatively narrow stress plateau. This is one of the un- 
satisfactory approximations to the GLAMM model that 
may explain why such a large Z value (compared to ex- 
periments) is required for the RP model to exhibit veloc- 
ity recoil and transient banding. This is an obvious and 
important direction for future work. 
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E. Constitutive Relations 

The RP model readily admits non-monotonic consti- 
tutive curves that can display shear banding, depending 
on the values of the CCR parameter j3, the entanglement 
number Z, and the solvent viscosity e. For smaller e and 
larger Z there is a wider separation between the high and 
low shear rate branches, while the CCR parameter con- 
trols whether the crossover between the two regimes con- 
tains a decreasing stress (smaller (3 or less active CCR) or 
is monotonic (larger (3 or more active CCR). Fig.[T]shows 
the effect of changing the number of entanglements for 
fixed solvent viscosity and CCR parameter, while Fig. [2] 
illustrates the effect of varying the the CCR parameter. 
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FIG. 2. Monotonic and non-monotonic constitutive curves 
arising from the Rolie-Poly model obtained by varying the 
CCR parameter /3, for e = 10~ 5 and Z = 265. The marked 
point at log 10 7T,j = 1.138 (7^ = 13.74) on the monotonic 
constitutive curve is studied in detail in Figs. (|4I5I6[) . 

The parameters space (/?, e, Z) for the Rolie-Poly 
model can be divided into regions of monotonic and non- 
monotonic constitutive curves, illustrated in Fig. [3] Non- 
monotonic curves occur for small values of e, small values 
of (3, and large Z. 

In what follows we will illustrate the transient banding- 
like behavior associated with monotonic constitutive re- 
lations. We will study the dynamics of startup to a 
shear rate that is on the plateau region, and which even- 
tually yields a non-banded homogeneous steady state. 
Transient banding behavior develops under these con- 
ditions, and we will explain this with the help of lin- 
ear stability analysis. We primarily use the parameters 
(/3,e, Z) = (0.71, 10 -5 , 265), corresponding to the thick 
line in Fig. [51 Unfortunately, Z — 265 is much higher 
than the values Z > 40 — 50 needed to see shear band- 
ing experimentally. We believe that this is because the 
Rolie-Poly model, as discussed above, only roughly ap- 
proximates the full microscopic theory of the GLAMM 
model, which is itself an approximation; hence we caution 
the reader against over-interpreting the values of these 
parameters. We will return later to briefly discuss this 
point. 
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FIG. 3. Lines separating monotonic (upper right) from non- 
monotonic (lower left) constitutive behavior as a function of 
the viscosity ratio e and CCR parameter /3, for different en- 
tanglement numbers Z. 



III. DYNAMICS AND BOUNDARY CONDITIONS 
A. Flow geometry 

We consider Couette flow between concentric cylinders 
of inner (outer) radii R\ {R2), in cylindrical coordinates 
with flow varying only in the radial direction, v = v(r)0. 
Using this description the tangential and radial com- 
ponents of the force balance conditions are respectively 
given by 



1 



pv = -jd r {r a r9 ) 



-p — = d r a rr 



(10) 
(11) 



and the components of the polymer strain W obey 
CW ee - 2W r6 i 

tr r- 

W rr ) 



— (l-A) + ^-(W rr ~Wgg) (12a) 



CWr: 



2 , „ 22?, 

(l-A) + —(Wg, 

tr r 2 



CW Z 



(1-/1) 

tr 

42? 

CW r e ~ Wrrl = 7 ^W r 9, 



(12b) 
(12c) 

(12d) 



where 7 = rd r (-) is the local shear rate and we have 
defined the nonlinear operator 



V 

C = dt ctrct 



1 

Td T R 



[I- A (W)] [1 + 13 A (W)] . 

(13) 

The constitutive equations can be simplified by 
parametrizing the spatial coordinate by 



r = R ie q 



(14) 



where 



6 



1 R * 
9 = ln i?T 



(15) 



is a measure of the curvature of the measurement de- 
vice, and the dimcnsionless coor dinate y € [0, 1 ] span s 
the width of the gap L = R 2 -R\ (|Greco and Ball 1 19971) . 

The relative stress difference between the two cylinders 
is given by 



oYe(-Ri) - cr rS (R 2 ) 



= 1 - e 



-2« 



(16) 



A Couette cell with radii R\ = 2 cm, R 2 = 2.1cm has 
a = 0.049. In a cone and plate device the stress varies as 
(lLarsod . fl"999h 



a g<t> {9 = tt/2) 



(17) 



where 9 is measured with respect to the normal vector 
to the plate. Hence the stress variation between the cone 
and plate is 



tan 2 a, 



(18) 



where a is the cone angle. For cone angles a = tt/2 — 9 = 
4° and 1° the stress difference is roughly equivalent to 
the values q = 2 x 1 0~ 3 and q = 2 x 10~ 4 , respectively 
(| Adams et a/J . l2008h . _n what follows we will sometimes 
compare the stress variation of a cone and plate geometry 
to the equivalent stress variation of a Couette cell with 
a particular geometric parameter q: however, we do not 
perform calculations for the cone and plate geometry. 

We change to dimcnsionless quantities, labelled by ~, 
by measuring time in units of the disengagement time Td 
and stress in units of the plateau modulus G, so that the 
total shear stress is expressed as 



O r e = W r 6 + £7, 



where 7 = jra and a r g 
less diffusion constant 



a r e/G. We define a dimension- 



D 



(20) 



Note that the inner radius R\ can be written in terms of 
the gap L between cylinders as 



Ri 



L 



e« - 1' 



(21) 



so that qR\ depends only weakly on q for small q. Below 
we will use D and q as independent parameters. Because 
we consider small q < 10~ 2 , the dependence of D on 
q is negligible. Hence, by fixing D and varying q we 
can separate the effects of the total stress gradient, as 
parametrized by q, from the effects of the finite width 
\fD of the interface. 



In this representation the constitutive equation be- 
comes 

CX = -6(1 - A)Z + 2^7 + 2De~ 2qy (Y - X) (22a) 
CY = -6(1 - A)Z - 2De~ 2qy (Y — X) (22b) 
CW =-&{!- A)Z (22c) 
CS= (1 + Y)tf -4De- 2qv S, (22d) 

where X = Wee, Y = W„, W = W zz , and S = W r e, 
and Z is the entanglement number given by Eq. ([5]). The 
nonlinear operator is thus 



C = d f 



De- 2q yd 2 + l + 6Z(l-A){l + (3A), (23) 



and the momentum balance equations become 



pv 
-pv 



e-^d y (e 2 ™a re ) 



In - 

Oy(J rr 

q 



+ (a rr - dee)- 



(24) 
(25) 



where p = f^. The dimensionless velocity is scaled 

by Td and a length R\(e q — 1) that becomes the plate 
separation R 2 — R\ = L in the q — > limit. 

Planar Couette flow obtains in the limit q — > 0, with 
the correspondences f — !> — y, 6 — > x, q — > and qR\ — > 
L, and the force balance conditions become pv = d y (a yx ) 

and dyOyy = 0. 



B. Boundary Conditions 

An imposed cylinder rotation rate V/Ri can be ex- 
pressed in terms of the integral of the shear rate across 
the gap by 



V 
R~i 



B 2 



dr 



Ri 



= / jqdy 



(19) 

The average applied shear rate is 



<7) = 



dr 



= idy. 



(26) 



(27) 



which reduces to —V±/L in the planar limit. We assume 
that there is no slip at the walls, despite the fact that 
wall slip can b e important experimentally (IWand . Il999t 
iBoukanv et~ai\ , 120061: iBoukanv and Wand . l2009af) . 

The presence of spatial gradients in the form of 
the 'diffusion' term necessitates a boundary condition 
on the viscoelastic strain W. Although it is possi- 
ble to incorporate sophisti cated wall constitutive models 
( Black and Graham . 19961 ) or complex bounda ry condi- 
tions (|Rossi et all l2006p|Adams et all 120081) . we will 
use the simplest Neumann boundary conditions for the 
polymer strain, 



(n ■ V)W = for y = 0, 1; 
here n is the normal to the boundary. 



(28) 



7 



IV. CALCULATION AL METHODS 



(Eq. Q3D is 



A. Linear Stability Analysis 

Wc first analyze the linear stability of Eqs. (|2"2"j) above, 
in the planar limit (q — > 0), which should be a good 
approximation of the behavior of typical systems for 
which q ~ 0.0001 - 0.05. Wc follow the stability 
analysis performed for theJohnson-Segalman model by 
iFielding and Olm sted (I2003D . who showed that fluids 
with non-monotonic constitutive curves are unstable dur- 
ing startup at fixed shear rate to spatial perturbations, 
which ultimately develop into shear bands. We will ap- 
ply this to monotonic constitutive curves such as shown 
in Fig. [U 

The equations of motion are first solved for a homoge- 
neous time-dependent solution, which we refer to as the 
homogeneous base state. We study the time dependence 
of inhomogeneous perturbations about this homogeneous 
base state. If these perturbations grow in time then the 
homogeneous time-dependent base state is unstable to in- 
homogeneous states, which we refer to as transient band- 
ing. Hence, we linearize the Navier-Stokes equation and 
the Rolie-Poly constitutive equations around a uniform 
time dependent base state (denoted by subscripts 0), 
X = _X + 6X, Y = Y + 5Y, W = W + SW, S = S + 5S 
and 7 = 70 + S-y. The linearized Navier-Stokes equation 



pdtdj 



(29) 



where 7 = d y v. Note that p is exceedin gly small (p ~ 
10~ 9 ) for the experiments of interest ( Tapadia et all 
120061 ) , and justifies the creeping flow limit in section llVBI 
However, we keep this term in the linear stability anal- 
ysis for completeness; leaving it out has no discernible 
effect for physically realistic values. 

We now consider spatial fluctuations in all quantities, 



iky 



u = (X,Y,W,S,% 



(30) 
(31) 



where Suk — <5uL fc because the function Su(y) is real. 
The boundary conditions of Eq. (|2"o) on the shear rate 
7 quantize the values of the wavenumber, k = mr (n = 
1,2,3,...), to keep the average shear rate across the gap 
fixed to the imposed value. 

Linearizing Eqs. (|22[ |2"9"|) about the homogeneous in- 
stantaneous state produces a matrix equation of the form 



d t Su k = M(fc,t)-*Ufc, 
where the stability matrix is given by 



(32) 



M 



n- zA 3 (i + rx) -A 3 z{i + rx) 



-a 3 tzw 

~A 3 TZS 




-a 3 tws 




-A 3 Z{\ + TX) 


21 


2S 


-A 3 Z(1 + TY) 








- ZA 3 (1 + YW) 








-A 3 TZS 


n 


1 + Y 





_ *d 
P 


p 



(33) 



with 



perturbed state at a given time. 



n = -1- Dk 2 -6Z(1- A)(l 
r = 1 - p + 2Ap. 



■Ap) 



(34) 
(35) 



The matrix depends on the values of the shear rate 
and polymer strain at a time t, determined according to 
homogeneous evolution of the dynamics. A linear distur- 
bance <5ufc will grow exponentially in time if an eigenvalue 
ma o f M has positive real part ( Fielding and Olmstedl . 
120031 ). We define w max as the eigenvalue with the largest 
real part 3?(a; Q ). To analyze the stability of the tran- 
sients with respect to spatial fluctuations we calculate 
the evolution of the eigenvalues of M as a function of 
time t, using the homogeneous solution W(i) for the un- 



A useful measure of the duration and severity of in- 
stability is the total 'weight' f2 max of the instability, ob- 
tained by integrating the real part of w max over all time 
during which it is positive: 



E 



Max(0, »(aw))tft. 



(36) 



For larger weights fi max we expect the fluid to enter the 
non-linear regime and develop a transient shear band, 
after which stable homogeneous flow is restored in steady 
state (for monotonic constitutive curves) . This is verified 
by the full non-linear calculation, to which we now turn. 
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B. Spatially resolved Rolie-Poly model 

We have thus far addressed the evolution of the homo- 
geneous model, and the linear (in)stability of this evo- 
lution to inhomogcncous states. To compute the full in- 
homogeneous non-linear dynamics we solve the creeping 
flow equation, Eq. and the DRP constitutive equa- 
tions, Eqs. (|22j). with a constrained shear rate, Eq. (|26|) 
and the boundary conditions, Eq. (|2"8|) . This leads to a 
set of coupled second order partial differential equations. 
A uniform spatial grid of 500 points was used to dis- 
cretize the equations, and a semi-implicit Crank-Nicolson 
scheme was used to effect the time evolution. We have 
checked convergence with respect to timcstep and spatial 
mesh. Typical dimensionless time steps were 10 -5 , with 
dimensionless spatial meshes of 2 x 10 -3 ; and a dimen- 
sionless diffusion constant D = 4 x 10~ 4 was used, cor- 
responding to an effective dimensionless diffusion length 
2 x 10~ 2 . 

The shear rate constraint of Eq. ([26]) was implemented 
by integrating Eq. (|24[) in the zero Reynolds number limit 
to obtain 



Te~ 2qy = a r6 = W, 



T 



1 - e-' 2 i 



[(W rB ) + e (7>] , 



(37a) 
(37b) 



where T is the torque per unit length applied to the inner 
cylinder. Eqs. (|3"?| were then used to eliminate the local 
shear rate j(y ) in Eqs. (|22[) in terms of the torque and 
the local value of the polymer shear strain W r e(y): 



1 



}(y) = - [Te- 2 ^ - Wre{y)) 



(38) 



The torque per unit length T can be computed at any 
instant in time from the full inhomogeneous profile in 
terms of the spatial averages of W r g and 7, as defined in 
Eq. (I37E1) . 

There are several possible measures of the degree of 
inhomogeneity in the spatially resolved flow profile. At 
any instant we define the shear rate drop 



7n 



7n 



(7) 



(39) 



where 7 max (rcsp. 7 m i n ) is the maximum (resp. mini- 
mum) of the measured shear rate profile, which should 
be easily accessible from experimental data. Other mea- 
sures could include the variance of the shear rate values 
across the sample cell, or the maximum local gradient of 
the shear rate. Since all of these give an equivalent qual- 
itative indication of the development of inhomogeneities 
during flow, we choose this simple measure. 



V. RESULTS: INSTABILITY AND TRANSIENTS 

In this section we study the linear (in)stability of 
homogeneous startup, and compare this with the non- 
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FIG. 4. Evolution of various quantities as a function of 
time during shear rate startup for parameters (/3, e, Z) = 
(0.71, 10" 5 ,265), for q = 10" 3 and applied average shear rate 
jTd — 13.74 (log 10 jTd — 1.138) (the marked point on the 
monotonic constitutive curve in Fig. [2]). (i) Dimensionless to- 
tal shear stress calculated based on homogeneous flow (solid) 
or allowing inhomogeneous transients (dashed); (ii) real part 
of the largest unstable eigenvalue t<j max ; (iii) shear rate drop 
A (as defined in Eq. I39[) , as a measure of spatial inhomogene- 
ity showing significant transients during startup. The times 
a, b, c, d, e are shown as velocity profiles in Fig. [5] 



linear spatially-resolved calculations in which the tran- 
sient shear bands are allowed to form naturally. Here 
we consider an initial state at rest with no 'noise', in a 
weakly curved geometry q — 10 -3 . In Sec. I VII we will 
study the effects of different magnitudes of spatial noise 
and curvature. 

We study model parameter values (/?, e, Z) = 
(0.71, 10 -5 , 265), corresponding to the solid constitutive 
curve in Fig. [2] for an imposed shear rate JTd = 13.74 
(the marked point at log 10 7Td = 1.138) and with shear 
stress gradient parametrized by q = 10~ 3 . During 
startup the shear stress overshoots, before settling down 
to steady state on a time of order tj,. Fig. [4] compares a 
calculation in which the shear rate remains homogeneous 
(solid line) with a spatially-resolved calculation that al- 
lows for inhomogeneities (dashed line), showing that 
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the stress decays more rapidly in the spatially-resolved 
model. The largest unstable eigenvalue, from the lin- 
ear stability analysis, becomes positive for log 10 (t/r c ;) > 
— 1.2, and reaches a maximum shortly after the stress 
overshoot. At late times it decays to a slightly negative 
value, controlled by the very shallow slope of the steady 
state constitutive curve. 
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FIG. 5. Velocity profiles during startup flow for the param- 
eters of the marked point log 10 jt<i = 1.138 (77^ = 13.74) 
on the monotonic constitutive curve in Fig. [5] (/?, e, Z) = 
(0.71, 10 _5 ,265), for q = 10~ 3 . The corresponding ellipsoids 
indicating the components of the W rr , Wee and W r e com- 
ponents of the polymer strain are shown below. At time 
t = 0.25r,j, which is shortly after the stress overshoot, both 
negative recoil velocities and transient shear banding occur. 
The times labeled a, b, c, d, e are the points labeled on A(t) 
in Fig. H 

Associated with this positive eigenvalue is a pro- 
nounced inhomogeneous transient during startup, as ev- 
idenced by the shear rate drop A in Fig. 2] and the veloc- 
ity profiles in Fig. [5] The most significant heterogeneity 
comes shortly after the eigenvalue has reached its maxi- 
mum, after the instability has had sufficient time to de- 
velop non-linear consequences. Fig. [5] shows that this 
maximum (at time t = 0.257^, labeled b) corresponds to 
a dramatic transient banding profile, in which the veloc- 
ity field has recoiled and become negative at the outer 
part of the cell (closer to y — 1), where the total stress 
is slightly lower. The strong transient bands decay af- 
ter a few Td until the steady state profile is reached. In 
this case, the stress gradient corresponding to q = 1CT 3 is 
enoug h to induce significant hetero geneity even in steady 
state ( Adams and Olmstedl . l20"09bl ) (we will return to the 
role of the stress gradient in Sec. IVI[) . Fig. [S]also shows 
the degree of polymer deformation (strain) during the re- 
coil. The polymer strain tensor becomes well-ordered, at 
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FIG. 6. The weight f2 max of the instability (magnitude of 
the integral of the most unstable eigenvalue) and the mono- 
tonic constitutive curve (Jxyi^f) (solid line) from Fig. [2] For 
03, e, Z) = (0.71, 10~ 5 , 265), for q = 10~ 3 . 



an angle of approximately 45° with respect to the flow 
direction, when the transient develops (b). Interestingly, 
it remains ordered and changes only weakly until steady 
state is reached; moreover, the spatial gradient only in- 
troduces a relatively small gradient in the molecular or- 
der and orientation. Hence, the transient inhomogenc- 
ity would be difficult to observe using a molecular probe 
such as birefringence or neutron scattering, as compared 
to explicit measurements of velocity profiles. 

Having explored the response to a single shear rate 
startup protocol, we now study the strength of the insta- 
bility for different shear rates. Fig. [5] shows the weight 
!] max of the instability (integrated unstable eigenvalue) 
for all applied shear rates for this parameter set. The 
largest weight occurs near the flattest part of the stress 
plateau, and indeed we find the clearest transient band- 
ing in this region. A detailed study (not shown here) 
reveals that the weight f2 max is larger for large Z, small 
e and larger /3; i.e. parameters for which the two flow 
branches are more widely separated and the constitutive 
curve has a shallower slope. The larger value of /3 pro- 
duces a flatter, longer stress plateau. 

In summary, transient banding is initiated by a lin- 
ear instability that begins during startup, and becomes 
largest shortly after the stress overshoot. If this in- 
stability has enough time and strength to grow, as 
parametrized by the weight O ma x, then transient band- 
ing and even recoil can appear before the eventual steady 
state is reached. The steady state is not banded, but 
weakly inhomogeneous as specified by the stress gradient 
across the rheometer. 



VI. THE ROLE OF PERTURBATIONS IN TRANSIENT 
BANDING 

We have seen thus far that inhomogeneous flow, or 
cquivalently transient shear banding, is correlated with 
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the most unstable eigenvalue governing fluctuations away 
from a homogeneous state at a given instant in time. We 
now study the spatial perturbations needed to trigger 
the instability, by evolving the non-linear inhomogenous 
equations of motion. There are several sources of spa- 
tial inhomogeneity that could serve this purpose: curved 
rheometer devices (cone and plate or cylindrical Couette 
geometries) possess a spatial inhomogeneity in the shear 
stress; thermal or instrument noise is presumably always 
present; thermal gradients across the cell may persist; 
and there may be inhomogeneous initial conditions due 
to very slow relaxation after sample loading. Here we 
will study (i) the consequences of a stress gradient due 
to curved streamlines and (ii) the role of inhomogeneous 
initial conditions that may arise due to the other effects 
just mentioned. 

We first consider the effect of random initial conditions 
in a geometry with no intrinsic stress gradient (i. e. pla- 
nar Couette flow, q = 0). A cosine wave with wavelength 
2L, in keeping with the boundary conditions, was set 
as the initial condition for a given polymer strain com- 
ponent, with the others set to zero. The constitutive 
equations were then evolved and the velocity profile cal- 
culated. This was carried out for a range of different am- 
plitudes A of the cosine wave. The different components 
of the polymer strain have markedly different effects: 
perturbations in W yy have the most dramatic effect, be- 
cause this strain component (along the flow gradient) is 
advected into the polymer shear strain W xy (Eq. I22dl) . 
which in turn directly contributes to the measured total 
shear stress. Conversely, perturbations in W xy are much 
less important because they are rotated into the velocity 
direction (W xx ), which does not contribute to the total 
shear stress. Fig. [7] shows the effect of an initial con- 
dition in W yy . For small enough initial noise there is 
virtually no transient banding, while substantial initial 
noise (a few percent of the steady state polymer strain 
in flow) can lead to strong transient banding and nega- 
tive velocity recoil. The modulus of entangled polymers 
is G ~ cksT, where c is the concentration of entan- 
glement strands. Reasoning that coherent fluctuations 
can obtain wit hin a volume of order a 3 , whe re a is the 
tube diameter, M arrucci and GrizzutH (]1983[) estimated 
the typical average strain (SWap) due to thermal fluctu- 
ations in entangled polybutadiene melts to be of order 
y^l/6.4 ~ 0.39, assuming a ~ 5.8 nm; this is comparable 
to the initial condition necessary to induce a transient 
banding response in our numerical calculation. 

The stress overshoot also changes character (Fig. [8]): 
for a larger value of A (a stronger perturbation) the stress 
relaxes more quickly after the overshoot, due to the more 
rapid development of transient banding. For A = 10~ 2 
the stress actually undershoots the steady state value be- 
fore increasing and then finally decreasing slowly to the 
ste ady state. Similar s tress undersh oots were reported 
by iTapadia and Wanel (|2004 ) and by ISui and McKennal 
(|2007l ). for entangled p olymer solutions that underw ent 
flow instability; and bv lCrawlev and Graesslevl ( 19771 ) in 



much earlier work. 

An alternative perturbation is a stress gradient due to, 
e.g., the curvature of a Couette cell or a non-zero angle 
in a cone and plate rheometer. This is studied in Fig. |H1 
As explained earlier we parametrize this by the curva- 
ture parameter q = I11R2/R1 (Eq. [T5l. For q > 10~ 3 , 
whose stress gradient is similar to that of very thin gap 
Couette rheometers or typical cone and plate rhcomctcrs, 
there is substantial transient banding for these parame- 
ters; while for q = 10 -5 there is essentially no transient 
banding. [We remind the reader that this calculation was 
performed for cylindrical Couette flow, and that the com- 
parison with cone and plate flow is based solely on the 
shear stress gradients in the two geometries.] The stress 
overshoot also relaxes quicker after the overshoot for a 
larger stress gradient (larger q) fFig. ITU]). 

An experimental test is thus to observe the change in 
transient banding as a function of flow geometry and 
hence q; or to systematically pre-shear the material to 
induce molecular deformation; or otherwise induce or 
spatially inhomogeneous initial conditions into a slid- 
ing plate geometry. Existing experiments show transient 
banding in both curved and flat geometries; our calcu- 
lations would be consistent with reproducible transients 
in a curved geometry, but less reproducible transients in 
a flat geometry, reflecting the variations of initial and 
loading conditions. 



VII. DISCUSSION 

A. Parameter Space and Critique of the Rolie-Poly Model 

We have explored the behavior of transient banding as 
a function of the parameters of the DRP model (Fig. [TT|l . 
A monotonic constitutive curve is more (dynamically) 
unstable and susceptible to transient banding if it has 
well separated flow branches, which occurs for greater 
numbers of entanglements Z, larger convective constraint 
release (CCR) efficiency j3, and lower solvent viscosity 
e. The larger CCR parameter j3 may appear paradox- 
ical, since increasing it renders a banding fluid stable. 
However, one also needs, simultaneously, a smaller vis- 
cosity ratio and larger Z, to render the plateau shal- 
lower and wider in the monotonic regime. Although 
the value of the CCR parameter /3 is evidently impor- 
tant, we have as yet no way of indep endently determin- 
ing it. iLikht man and Grahaml (|2003t ) used the values of 
,3 = 0, 0.5 and 1 to fit multimode, steady state and tran- 
sient data respectively. It would be useful to compare 
a quantitative model such as the DRP model with ex- 
perimental results on shear banding (transient or steady 
state) to try and infer the magnitude of the CCR param- 
eter, given knowledge of the other parameters. 

Unfortunately the DRP model is still too crude to 
make this a justifiable exercise. The high shear rate 
branch is not handled correctly, since Newtonian behav- 
ior is not observed experimentally at the highest shear 
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FIG. 7. Effect of noisy initial conditions for W yy on (a) the shear rate drop A as a measure of inhomogeneity, and (b-d) 
transient velocity profiles v(y), for different initial conditions W yy = Acos7ri/. In steady state W yy = —0.52. For parameters 
(0, e, Z) = (0.71, 10~ 5 , 265) and q = 0. 
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Acos7ry, parameterized by A. For parameters 



(P,e,Z) = (0.71, 10~ 5 ,265). 



rates (|Tapadia and WangL 120031 ). Moreover, the DRP 
model is only a simple a pproximation to the mor e com- 
plete GLAMM the ory (jLikhtman and Grahaml . 120031: 
iGraham et aill2003l) . and there is thus not a precise cor- 
respondence between the parameter values used to fit 
the DRP model to data, and their microscopic physical 
meaning in the original GLAMM model. For example, we 



use rather larger values of entanglement number Z than 
are realized in most experiments. However, the DRP 
model reproduces the qualitative behavior, albeit with- 
out the quantitative material parameters. These quali- 
tative conclusions seem to be robust, as far as we can 
discern: the much broader and flatter steady state con- 
stitutive curves are accompanied by dramatic dynamics 
that can result in transient inhomogeneitics. There is 
an urgent need for performing inhomogeneous calcula- 
tions with more detailed theories such as the GLAMM 
model. The GLAMM model itself is also only an ap- 
proximation, and it certainly misses important physics 
about the behaviour of the tu be at very high shear rates 
(|Graham and McLeishl . l2007t ). 

Fig. [TT] shows the behavior expected for three values 
of /3, as a function of entanglement number Z and vis- 
cosity ratio e, for q = 10 -3 . The figure shows regions 
of (1) steady state banding signified by a non-monotonic 
constitutive curve; (2) transient banding during startup; 
and (3) dramatic recoil with negative velocities during 
startup. From our startup calculations we found these 
phenomena to occur for shear rate drops (see Eq. 131?)) 
A > 2 (transient banding) and A > 7 (transient banding 
featuring srecoil). Experiments show transient banding 
for Z ~ 15 — 50. steady state banding for higher Z, and 
recoil for some fluids that do not shear band in steady 
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FIG. 9. Effect of different stress gradients q = In -g* on (a) the shear rate drop A as a measure of inhomogeneity, and (b-d) 
transient velocity profiles v(y). For parameters (/3, e, Z) = (0.71, 10~ 5 , 265). The value q = 10~ 3 corresponds roughly to a cone 
angle of 2.5°. 
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state. Based on this, we suggest that the experimental 
phenomenology of entangled polymer solutions most re- 
sembles that for values of the CCR parameter (3 in the 
range between /3 = 0.71 and /3 = 0.8, as parametrized by 
the Rolie-Poly model. We emphasize that the parameter 
values for this model do not necessarily correspond to 
those of experiments. However, polymer solutions have a 



reasonably wide range of Z for which transient, and not 
steady state, banding is found as exhibited in Fig. [TT] for 
f3 = 0.71. 

B. The stress overshoot and the non-monotonicity of the 
instantaneous constitutive curve 

Fig. @] shows that transient banding occurs during the 
decreasing stress after the stress overshoot, and occurs 
a short time after the unstable eigenvalue is most un- 
stable. At these shear rates the fluid behaves elastically 
at early times, and the stress overshoot as a function of 
time could thus be envisioned as the stress overshoot of 
a solid as a function of strain 7 = jt] such an elastic ma- 
teri al with da xv /d^ < could be expected to be unsta- 
ble jMar^cci^dGri^tl Il983t) . This was also noted 
bv lSui and McKennal (|2007h in their study of entangled 
polymer solutions 

Fig. [S] shows that the integrated instability, as deter- 
mined by the linear stability analysis, is most pronounced 
near the flattest part of the constitutive curve. Since the 
steady state is nearly unstable at this point, it is tempting 
to declare the nearly flat slope as the cause of the insta- 
bility. However, the steady state behavior may have lit- 
tle influence on the instantaneous dynamics. To explore 
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The numbers in boxes denote the largest value of the shear rate drop A found during startup for an applied average shear 
rate on the flattest part of the steady state constitutive curve, analogous to the solid circle in Fig. [6] Shown are: steady state 
banding in the non-monotonic region; transient banding (A > 2 and diamonds or ellipses); and negative velocity recoil upon 
startup (A > 7 and ellipses). 



this, we have calculated the instantaneous constitutive 
cu rve, as descr i bed in the introduction and constructed 
bv lHaves et all ( 2010f ) from their experimental data. We 
calculate the startup flow for different homogeneous shear 
rates to obtain a surface a xy (j,t). This surface then de- 
fines a curve relating shear stress to shear rate at any 
instant in time t (the instantaneous constitutive curve). 
If this curve has a negative slope, then one expects an 
instability in the dynamics at finite Reynolds number, 
or at least in a subspace of the full dynamics specified 
by the matrix M that appears in the linear stability 
calculation of Section II V A I FigJT^] shows that the in- 
stantaneous constitutive curve becomes non-monotonic, 
with the strongest instability (largest unstable eigenvalue 
w m ax) indeed occurring very close the point of steepest 
negative slope daxy/d'y. 

Hence, transient banding, while evidently linked to the 
steady state constitutive curve, may be expected when 
the instantaneous constitutive curve is unstable (has neg- 
ative slope), based on linear stability analysis. It may 
also be expected during a stress overshoot, for applied 



shear rates that approximate a step strain and are thus 
strong enough to bring the fl uid into its elastic regime 
( Marrucci and Grizzutil . Il983l ). We emphasize that we 
have shown a very strong correlation, but in the context 
of the Rolie-Poly model. The precise correspondences re- 
main to be explored in other models, and we await more 
work on this subject. 



C. Other transient phenomena 

Two other transient banding phenomena have been 
reported. Well defined shear bands have been ob- 
served during portions of the cycle during Large Am- 
plitude Osdlaiory Shear ( LAOS) in polymer solu- 
tions ( Tapadia et all 120061 ; iRavindranath and Wand . 
2008). This has al s o been calculated numerically by 
Adams and Olmsted! (|2009bf ) using the DRP model for a 
monotonic constitutive curve, and by IZhou et~cd\ (120101 ) 
for a non-monotonic constitutive model devised for 
wormlike micelles. The calculations revealed banding- 
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FIG. 12. (a) Maximum eigenvalue as a function of shear rate 
7 for different times t after startup, (b) instantaneous consti- 
tutive curves, calculated from starting up at different shear 
rates 7. Each curve represents the stress obtained a time t 
after starting up at the given shear rate 7. Note the strong 
correlation between the region of positive maximum eigen- 
value, and the region of negative slope of the time dependent 
constitutive curve. 



like transients for shear rates sweeping across the pseudo- 
plateau of a barely monotonic constitutive model, and for 
frequencies slower than the inverse reptation time. This 
is slow enough for a transient band to form during the 
upward sweep of strain rate (according to Fig. [4]), but 
not so slow that it decays away completely. It can then 
be reinforced in subsequent sweeps, given the strong spa- 
tial perturbation it represents, until it develops a sharp 
profile after a few cycles. 

In other experiments, step strain at shear rates 777? > 
1 leads to an inhomogeneous response (including nega- 
tive veloci ty recoil) aft e r cess a tion of step strain, in both 
solut ions ( Wang et ai . 20061 IRavindranath and Wand . 
120071 ) and melts ( Boukanv et aLl . l2009l ). In one example, 
a solution of Z = 64 was strained to just the stress max- 
imum in the overshoot. During relaxation, recoil and an 
inhomogeneous r esponse developed, similar to our previ- 
ous calculations ( Adams and Olmsted . 2009a). One ex- 
planation for this result is that the fluid became unstable 
during strain, and the instability was able to amplify a 
nascent fluctuation large enough to induce a non-linear 
and heterogeneous response during the subsequent re- 
laxation. Indeed, the calculation in Fig. |4] shows that 
the largest eigenvalue can become positive, signifying in- 
stability, before the overshoot is reached. However, this 
cannot be the entire story. In melts the inhomogeneous 
respo nse can also be a dramatic fracture (jBoukanv et 



could be because the DRP model inadequately captures 
the physics of the GLAMM model; or because impor- 
tant physics governing the heterogeneity associated with 
retracti on has yet t o be correctly incorporated in any 
model (|Wang et all l2007t) . 



D. Implications for validating constitutive models 

A significant amount of literature, both experi- 
mental and theoretical, has been devoted to test- 
ing constitutive models against the initial stress tran- 
sients, including the overshoot and subsequent rela xation 
(jLikhtman and Graham! 120031: IWang et aZ.1 . 120061) . This 



is typically done by comparing the stress transient for 
shear rate startup to a calculation that assumes homoge- 
neous shear flow in a planar geometry. As is evident in 
Fig. |U the relaxation after the overshoot is faster when 
the fluid is allowed to become inhomogeneous, as physi- 
cally occurs during transient banding. Cone and plate ge- 
ometries are often used because unlimited strains can be 
applied and the weak stress gradient is usually hoped to 
be negligible. Unfortunately, the transient banding and 
hence its associated stress response can be pronounced 
even for very small stress varations, as parametrized by 
q. A planar geometry should exhibit less transient band- 
ing if there are no significant perturbations (Fig. [5]). We 
are not yet able to quantify this level of perturbation. 

Another common benchmark is the Doi-Edwards 
damping function h(j,t) = cr xy ('-f ,t) / 'cr xy (0,t), defined as 
the stress relaxation cr 3; j,(7,i) after a finite step strain 
7, normalized by the relaxation ^^(O,^) that would 
occur in the linear limit of zero strain. In strongly 
entangled syste ms anomalou s behaviour in hH , t) has 
been reported ( Osakil . Il993t IVenerusl . 120051) . corre- 
sponding to faster relaxation than predicted by DE 
theory As noted above, dramatic heterogeneities can 
also occur durin g relaxation after a step s train, in 
the DRP model (jAdams and Olmstedl. l2009bl) and ex- 
periments on entangled polym e rs (Wang et all [20061 : 
IRavindranath and Wangl . 120071: iBoukanv et all l2009h . 
These heterogeneities will have a signature in the stress 
response to a step strain. Thus, it is conceivable that 
transient heterogeneities could ac count for the anomalous 
damp ing function measurements iMarrucci and Grizzutil 
(|1983D . 



I2009D t hat, a s yet, we cannot calculate from the DRP 



model (jWand . 120091 ; I Adams and Olmstedl . [2009al) . This 



Hence, one must take great care in making predictions 
for constitutive models: perturbations such as noise or 
the stress gradients of curved geometries, or the presence 
of non-monotonic constitutive behavior, necessitates a 
validation against a full i nhomogeneous calc ulation. This 
warning was also given bv lZhou et al. I (|2008T) . Such calcu- 
lations are currently out of reach for models such as the 
GLAMM model, which already involve solving partial 
differential equations in time and arc length coordinates 
dGraham et all 120031) . 
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VIII. SUMMARY 

Experiments on highly entangled polymers dissolved 
in their oligomers show inhomogeneous transient velocity 
profiles in shear flow in a variety of experimental tests. 
We have modeled this behavior using the single mode 
diffusive Rolie-Poly (DRP) model. The DRP model has 
some quantitative shortcomings, but can explain many 
of the experimental transient phenomena. By varying 
the model parameters, namely the convective constraint 
release efficiency /3, the entanglement number, Z and 
the solvent viscosity e, we have shown that transient 
shear banding can occur even for mon otonic const i tutive 
curves. Similar behavior was found bv lZhou et all (|2008l ) 
in calculations using a different non-monotonic model 
for entangled polymers, the partially extended convec- 
tive strain model. 

A linear stability analysis of the startup transient 
shows that the homogeneous state is unstable, at early 
times, to the inhomogeneous transient shear bands. The 
"weight" f2 max of this instability, given by the time in- 
tegral of the most unstable eigenvalue, provides a good 
predictor of flow instability. Sl max is largest where the 
constitutive curve is flat. If this weight is large then 
small perturbations can provide the seed for this insta- 
bility. For large enough perturbations, particularly in 
fluctuations of the normal stress, the flow profile devel- 
ops a shear banding state, which then decays over a few 
relaxation times once the eigenvalue returns to a sta- 
ble value. We have shown that spatial perturbations play 
a crucial role in triggering the transient instability. Ex- 
amples include strong enough stress gradients, as found 
in rotational rheometers; or inhomogeneous initial condi- 
tions due to thermal noise; or residual stress from loading 
the sample. 

Hence, shear banding ranges from the steady state 
shear banding seen in non-monotonic constitutive curves, 
to transiently inhomogeneous flow that can develop even 
for monotonic constitutive curves. The stress overshoot 
during startup relaxes more quickly, due to the transient 
shear banding, than if the fluid were to remain homoge- 
neous. This may be consistent with the suggestion long 
ago by iMarrucci and Grizzutil ( 19831 ) that an elastic in- 
stability could lead to inhomgeneities in a step strain ex- 
periment. The relaxation of the stress overshoot is more 
pronounced in geometries with stronger stress gradients. 

Our calculations suggest that fluids with a non- 
monotonic instantaneous constitutive curve a xy (j,t) are 
more likely t o have an instability and transient banding 
( Haves et all l2010l ). Indeed, the linear stability anal- 
ysis in such situations leads to a unstable eigenvalue 
that is coincident (in time) with the more strongly non- 
monotonic instantaneous constitutive curves. Similarly, 
transient banding is strongly correla ted with the stress 
overs hoot observed during startup ( Sui and McKennal . 
I2007D . This is a subject for future work. 

In attempting to fit the parameters of the Rolie-Poly 
model to the existing data, it is apparent that the model 



should be judged mainly for its qualitative conclusions: 
for example, shear banding in the model occurs for much 
larger Z than is found in experiments. The model also 
fails to describe the high shear rate branch correctly with- 
out violating certain physical criteria used to define the 
effect of stretch on conv ective constraint release in the 
parent GLAMM mo del (|Likhtman and Grahaml . l2003t 
iGraham et cJll2003h . Finally, challen ges for future work 
include explaining the dramatic fracture seen in some re- 
cent experiments on b o th entangled melts and solution s 
([Boukany et all . 120091 : iRavindranath and Wanl 120071) . 
and extending (or replacing) the GLAMM model to ac- 
curately describe the highest shear rate behavior. 



ACKNOWLEDGMENTS 

We arc grateful to the UK EPSRC (SMF, 
EP/E5336X/2) and the Royal Commission of 1851 
(JMA) for financial support; and to Ron Larson for 
useful comments. 



Adams, J. M., Fielding, S. M., and Olmsted, P. D., "The interplay 
between boundary conditions and flow geometries in shear band- 
ing: hysteresis, band configurations, and surface transitions," J. 
Non-Newt. Fl. Mcch 151, 10f-118 (2008). 

Adams, J. M. and Olmsted, P. D., "Adams and Olmsted reply:," 
Phys. Rev. Lett. 103, 067801 (2009a). 

Adams, J. M. and Olmsted, P. D., "Nonmonotonic models are not 
necessary to obtain shear banding phenomena in entangled poly- 
mer solutions," Phys. Rev. Lett. 102, 219802 (2009b). 

Berret, J.-F., "Rheology of wormlike micelles: equilibrium prop- 
erties and shear banding transition," in Molecular Gels, edited 
by R. G. Weiss and P. Terech (Springer, Dordrecht, 2005) pp. 
235-275. 

Black, W. B. and Graham, M. D., "Wall-slip and polymer-melt 
flow instability," Phys. Rev. Lett. 77, 956-959 (1996). 

Boukany, P. E., Hu, Y. T., and Wang, S.-Q., "Observations of wall 
slip and shear banding in entangled DNA solutions," Macro- 
molecules 41, 2644-2650 (2008). 

Boukany, P. E., Tapadia, P., and Wang, S. Q., "Interfacial stick- 
slip transition in simple shear of entangled melts," J. Rheol. 50, 
641-654 (2006). 

Boukany, P. E. and Wang, S-.Q., "A correlation between velocity 
profile and molecular weight distribution in sheared entangled 
polymer solutions," J. Rheol. 51, 217 (2007). 

Boukany, P. E. and Wang, S.-Q., "Exploring the transition from 
wall slip to bulk shearing banding in well-entangled DNA solu- 
tions," Soft Matter 5, 780-789 (2009a). 

Boukany, P. E. and Wang, S.-Q., "Shear banding or not in entan- 
gled DNA solutions depending on the level of entanglement," J. 
Rheol. 53, 73 (2009b). 

Boukany, P. E., Wang, S.-Q., and Wang, X., "Step shear of entan- 
gled linear polymer melts: New experimental evidence for elastic 
yielding," Macromolecules 42, 6261-6269 (2009). 

Cates, M. E., "Nonlinear viscoelasticity of wormlike micelles (and 
other reversibly breakable polymers)," J. Phys. Chem. 94, 371— 
375 (1990). 

Cates, M. E. and Fielding, S. M., "Rheology of giant micelles," 

Adv. Physics 55, 799-879 (2006). 
Crawley, R. L. and Graessley, W. W., "Geometry effects on stress 

transient data obtained by cone and plate flow," J. Rheology 21, 

19-49 (1977). 

Denn, M. M., "Extrusion instabilities and wall slip," Ann. Rev. Fl. 
Mech. 33, 265-287 (2001). 



16 



Divoux, T., Tamarii, D., Barentin, C, and Manneville, S., "Tran- 
sient shear banding in a simple yield stress fluid," Phys. Rev. 
Lett. 104, 208301 (2010). 

Doi, M. and Edwards, S. F., The Theory of Polymer Dynamics 
(Clarendon, Oxford, 1989). 

Doi, M. and Edwards, SF, "Dynamics of concentrated polymer 
systems. Part 4.-Rheological properties," J. Chem. Soc, Far. 
Trans. 2 75, 38-54 (1979). 

Fielding, S. M. and Olmsted, P. D., "Kinetics of the shear banding 
instability in startup flows," Phys. Rev. E 68, 036313 (2003). 

Fukuda, M., Osaki, K., and Kurata, M., "Nonlinear viscoelasticity 
of polystyrene solutions. I. Strain-dependent relaxation modu- 
lus," J. Polym. Sci. Polym. Phys. 13, 1563-1576 (1975). 

Graham, R. S., Likhtman, A. E., McLeish, T. C. B., and Milner, 
S. T., "Microscopic theory of linear, entangled polymer chains 
under rapid deformation including chain stretch and convective 
constraint release," J. Rheol. 47, 1171-1200 (2003). 

Graham, R. S. and McLeish, T. C. B., "Emerging applications for 
models of molecular rheology," J. Non-Newt. Fl. Mech. 150, 11— 
18 (2007). 

Greco, F. and Ball, R. C, "Shear-band formation in a non- 
Newtonian fluid model with a constitutive instability," J. Non- 
Newt. Fl. Mech 69, 195-206 (1997). 

Hayes, K. A., Buckley, M. R., Qi, H., Cohen, I., and Archer, L. A., 
"Constitutive curve and velocity profile in entangled polymers 
during start-up of steady shear flow," Macromolcculcs 43, 4412- 
4417 (2010). 

Hu, Y. T., Palla, C, and Lips, A., "Comparison between shear 
banding and shear thinning in entangled micellar solutions," J. 
Rheol. 52, 379-400 (2008a). 

Hu, Y. T., Palla, C, and Lips, A., "Role of Electrostatic Interac- 
tions in Shear Banding of Entangled DNA Solutions," Macro- 
molecules 41, 6618-6620 (2008b). 

Hu, Y. T., Wilen, L., Philips, A., and Lips, A., "Is the constitutive 
relation for entangled polymers monotonic?." J. Rheol. 51, 275 
(2007). 

Islam, M. T. and Archer, L. A., "Nonlinear rheology of highly en- 
tangled polymer solutions in start-up and steady shear flow," J. 
Polym., Sci. Part B 39, 2275-2289 (2001). 

Jagla, E. A., "Towards a modeling of the time dependence of con- 
tact area between solid bodies," J. Stat. Mech. 2010, P06006 
(2010). 

Kolkka, R. W., Malkus, D. S., Hansen, M. G., and Ierley, G. R., 
"Spurt phenomena of the Johnson-Scgalman fluid and related 
models," J Non-Newt. Fl. Mech. 29, 303-335 (1988). 

Larson, R. G., "A constitutive equation for polymer melts based 
on partially extending strand convection," J. Rheol. 28, 545-571 
(1984). 

Larson, R. G., The Structure and Rheology of Complex Fluids (Ox- 
ford University Press, New York, 1999). 

Larson, R. G., Sridhar, T., Leal, L. G., McKinley, G. H., Likht- 
man, A. E., and McLeish, T. C. B., "Definitions of entanglement 
spacing and time constants in the tube model," J. Rheol. 47, 809 
(2003). 

Likhtman, A. E. and Graham, R. S., "Simple constitutive equation 
for linear polymer melts derived from molecular theory: Rolic- 
poly equation," J. Non-Newt. Fl. Mech. 114, 1 (2003). 

Likhtman, A. E., Milner, S. T., and McLeish, T. C. B., "Micro- 
scopic theory for the fast flow of polymer melts," Phys. Rev. 
Lett. 85, 4550-4553 (2000). 

Lu, C-Y. D., Olmsted, P. D., and Ball, R. C, "The effect of non- 
local stress on the determination of shear banding flow," Phys. 
Rev. Lett. 84, 642-645 (2000). 

Manning, M. L., Daub, E. G., Langcr, J. S., and Carlson, J. M., 
"Rate-dependent shear bands in a shear-transformation-zone 
model of amorphous solids," Phys. Rev. E 79, 016110 (2009). 

Manning, M. L., Langcr, J. S., and Carlson, J. M., "Strain localiza- 
tion in a shear transformation zone model for amorphous solids," 
Phys. Rev. E 76, 056106 (2007). 



Marrucci, G., "Dynamics of entanglements, a nonlinear model con- 
sistent with cox-merz rule," J. Non-Newtonian Fluid Mech 62, 
279 (1996). 

Marrucci, G. and Grizzuti, N., "The free energy function of the 
Doi-Edwards theory: analysis of the instabilities in stress relax- 
ation," J. Rheol. 27, 433 (1983). 

McLeish, T. C. B. and Ball, R. C, "A molecular approach to the 
spurt effect in polymer melt flow," J. Poly. Sci. B-Poly. Phys. 
24, 1735-1745 (1986). 

Mead, D. W., Larson, R. G., and Doi, M., "A molecular theory 
for fast flows of entangled polymers," Macromolccules 31, 7895- 
7914 (1998). 

Menezes, E. V. and Graessley, W. W., "Nonlinear rheological be- 
havior of polymer systems for several shear-flow histories," J. 
Polym. Sci.: Polym. Phys. 20, 1817-1833 (1982). 

Milner, S. T., McLeish, T. C. B., and Likhtman, A. E., "Micro- 
scopic theory of convective constraint release," J. Rheol. 45, 
539-563 (2001). 

Moorcroft, R. L., Cates, M. E., and Fielding, S. M., "Age- 
dependent transient shear banding in soft glasses," Phys. Rev. 
Lett. 106, 055502 (2011). 

Morrison, F. A. and Larson, R. G., "A study of shear- 
stress relaxation anomalies in binary-mixtures of monodispcrsc 
polystyrenes," J. Polym. Sci. Part B-Polym. Phys. 30, 943-950 

(1992) . 

Olmsted, P. D., "Perspectives on shear banding in complex fluids," 
Rheo. Acta 47, 283-300 (2008). 

Olmsted, P. D., Radulescu, O., and Lu, C.-Y. D., "The Johnson- 
Segalman model with a diffusion term in cylindrical Couctte 
flow," J. Rheology 44, 257-275 (2000). 

Osaki, K., "On the damping function of shear relaxation modulus 
for entangled polymers," Rheol. Acta 32, 429-437 (1993). 

Osaki, K. and Kurata, M., "Experimental appraisal of the Doi- 
Edwards theory for polymer rheology based on the data for 
polystyrene solutions," Macromolccules 13, 671—676 (1980). 

Ravindranath, S. and Wang, S.-Q., "What are the origins of 
stress relaxation in step shear of entangled polymer solutions?." 
Macromolecules 40, 8031-8029 (2007). 

Ravindranath, S. and Wang, S.-Q., "Large amplitude oscillatory 
shear behavior of entangled polymer solutions: Particle tracking 
velocimetric investigation," J. Rheol. 52, 341 (2008). 

Ravindranath, S., Wang, S.-Q., Olechnowicz, M., and Quirk, R. P., 
"Banding in simple steady shear of entangled polymer solu- 
tions," Macromolecules 41, 2663-2670 (2008). 

Rehage, H. and Hoffmann, H., "Viscoelastic surfactant solutions: 
model systems for rheological research," Mol. Phys. 74, 933 
(1991). 

Rossi, L. F., McKinley, G., and Cook, L. P., "Slippage and migra- 
tion in taylor-couette flow of a model for dilute wormlike micellar 
solutions," J. Non-Newt. Fl. Mech. 136, 79 (2006). 

Spenley, N. A., Cates, M. E., and McLeish, T. C. B., "Nonlinear 
rheology of wormlike micelles," Phys. Rev. Lett. 71, 939-942 

(1993) . 

Sui, C. and McKenna, G. B., "Instability of entangled polymers in 
cone and plate rheometry," Rheol. Acta 46, 877-888 (2007). 

Tapadia, P., Ravindranath, S., and Wang, S.-Q., "Banding in en- 
tangled polymer fluids under oscillatory shearing," Phys. Rev. 
Lett. 96, 196001 (2006). 

Tapadia, P. and Wang, S.-Q., "Yieldlike constitutive transition in 
shear flow of entangled polymeric fluids," Phys. Rev. Lett. 91, 
198301 (2003). 

Tapadia, P. and Wang, S.-Q., "Nonlinear flow behavior of entan- 
gled polymer solutions: Yieldlike entanglement-disentanglement 
transition," Macromolccules 37, 9083-9095 (2004). 

Tapadia, P. and Wang, S.-Q., "Direct visualization of continuous 
simple shear in non-newtonian polymeric fluids," Phys. Rev. 
Lett. 96, 016001 (2006). 

Venerus, D.C., "A critical evaluation of step strain flows of entan- 
gled linear polymer liquids," Journal of Rheology 49, 277-295 
(2005). 



17 



Venerus, D.C. and Nair, R., "Stress relaxation dynamics of an en- 
tangled polystyrene solution following step strain flow," Journal 
of rhcology 50, 59-75 (2006). 

Vinogradov, G. V., "Critical regimes of deformation of liquid poly- 
meric systems," Rheol. Acta 12, 273 (1973). 

Vrentas, C. M. and Graessley, W. W., "Study of shear stress re- 
laxation in well-characterized polymer liquids," J. Rheol. 26, 
359-371 (1982). 

Wang, S. Q., "Molecular transitions and dynamics at polymer/wall 
interfaces: Origins of flow instabilities and wall slip," in Poly- 
mers in Confined Environments, Adv. Poly. Sci., Vol. 138 
(Springer, Berlin, 1999) pp. 227-275. 

Wang, S.-Q., "Comment on "nonmonotonic models are not neces- 
sary to obtain shear banding phenomena in entangled polymer 
solutions"," Phys. Rev. Lett. 103, 219801 (2009). 

Wang, S-.Q., Ravindranath, S., Boukany, P., Olechnowicz, M., 
Quirk, R., Halasa, A., and Mays, J., "Nonquicsccnt relaxation 



in entangled polymer liquids after step shear," Phys. Rev. Lett. 
97, 187801 (2006). 

Wang, S-.Q., Ravindranath, S., Wang, Y., and Boukany, P., "New 
theoretical consideration in polymer rheology: Elastic break- 
down of chain entanglement network liquid state," J. Chem. 
Phys. 127, 064903 (2007). 

Zhou, L., Cook, L. P., and McKinley, G. H, "Probing shear-banding 
transitions of the vcm model for entangled wormlikc micelle 
model for entangled wormlikc miccllar solutions using large am- 
plitude oscillatory shearing (LAOS) deformation," J. Non-Newt. 
PI. Mech. 165, 1462-1472 (2010). 

Zhou, L., Vasquez, P. A., Cook, L. P., and McKinley, G. A., "Mod- 
eling the inhomogeneous response and formation of shear bands 
in steady and transient flows of entangled liquids," J. Rheol. 52, 
591-623 (2008). 



